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Abstract. We prove a generalization of Li-Yau estimate for a 
family of second order linear parabolic equations. We show that 
this generalization is sharp in the sense that the equality case is 
achieved by a family of fundamental solutions different from the 
Gaussian heat kernels. 

1. Introduction 



The Li-Yau estimate [TOj is the following sharp inequality of any 
positive solution p t of the heat equation p t = Ap t on a manifold with 
non-negative Ricci curvature: 

n 

(1-1) Alogp t >--. 

There are vast generalizations of the above important inequality to 
various other geometric evolution equations. For a detail account of all 
these generalizations as well as the importance of the inequality (II. ip . 
see [13]. 

There are generalizations of the inequality (II. ip to other second order 
linear parabolic equations under the, so called, curvature-dimension 
conditions (see for instance [3]). They are estimates of the form 

Tl 

(1-2) Llogp t >--, 

where L is a linear differential operator without constant term and pt 
is a solution of the equation p t = Lp t . 

However, one of the important features of the inequality (II .ip is its 
sharpness in the sense that equality of (II. ip is achieved by the following 
solution of the heat equation on the Euclidean space: 



1 ( \x 

Pt\ X ) = /,_,W9 eX P 



(47Tt) n / 2 * V 4t 

On the other hand, the following is a solution 

/ 27r(exp(2tfc)-l) V" /2 ( -k\x\" 
X3) Pt{x) = { kexp&k) ) 6XP Uexp(2tA;) 
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of the equation 

(1.4) p t = Ap t - k (x, Vp t ) = Ap t - ( V ( ||x| 2 ) , Vp t 



where k > is a constant. 

The solutions (11 .3p never satisfy the equality case of (II. 2p . Motivated 
by this observation, we prove the following generalization of (II. ip . For 
simplicity, we assume, through out this paper, that the manifolds are 
compact without boundary. This assumption can be removed by a 
localization procedure as in [TD1 ITTj . 

Theorem 1.1. Assume that 

(1) Rc > (Rc denotes the Ricci curvature of the underlying Rie- 
mannian manifold), 

(2) A(AU+ l\VU\ 2 ) < nk 2 

for some non-negative constant k. Then any positive solution pt of the 
equation p t = Ap t + (V£7, Vpt) satisfy the following estimate: 

(1.5) Alogp t + ~AU > -^coth(H). 

By letting U = and k goes to 0, we recover the estimate (11. ip . 
Note also that the solution (jl.3p achieves the equality case of ( II. 5p and 
the assumptions of Theorem 11.11 with U = — ||x| 2 . 

Li-Yau type estimates are also called differential Harnack inequalities 
since one can recover the Harnack inequality [12] by integrating the 
estimate (II. ip along geodesies. An analogue of this fact also holds true 
in our setting. However, instead of integrating along geodesies, the 
correct paths in this case are the minimizers of the following functional: 

(1.6) c s , t (x,y) = inf f \\^r)\ 2 + V(^r))dr, 

where the infimum is taken over all paths 7 : [s,t] — > M joining x 
and y. The idea of considering functionals of the form (I1.6P already 
appeared in [TO] . 



Theorem 1.2. Assume that 

(1) Rc > 0, 

(2) A(AU+ l\VU\ 2 ) < nk 2 

for some non-negative constant k. Then any positive solution pt of the 
equation p t = Ap t + (VU, Vpt) satisfies the following estimate: 

Pt (y) £ (*f)V ! «pf-jM«,») + !%)-W] 



Ps{X) 
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By letting U = and k goes to 0, we recover the following Harnack 
estimate. 



Corollary 1.3. [T21 [TO] Assume that Rc > 0. Then any positive solu- 
tion p t of the equation p t = Ap t satisfy the following estimate: 



p s [X) \S / 

It is also known that Corollary 11.31 recovers the heat kernel compar- 
ison theorem of Cheeger-Yau [S] if we let p t be the heat kernel and 
letting s goes to 0. The same principle also works for Theorem II . 21 

Theorem 1.4. Assume that 

(1) Rc > 0, 

(2) A{AU+ ±|W| 2 ) < nk 2 

for some non-negative constant k > 0. Then the fundamental solution 
p t (x,y) of the equation p t = Ap t + (VU, V pt) satisfies the following 
estimate: 
(1.7) 

MX ' V) ~ GvrsinMH) ) 2 ^ (~4 M ^ V) + U{v) " U{x)) 
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In the case U = — 2> ^ e cos ^ function is given by 

. . k\y\ 2 coth(H) 
co,t(0,y) = g «nt 

(see the proof of Theorem 13.51) and right hand side of (11. 7p becomes 
the fundamental solution (II. 3p . Therefore, all inequalities in Theorem 
11.41 become equalities in this case. 

Again, by setting U = and letting k goes to 0, we recover the 
Cheeger-Yau estimate. 

Corollary 1.5. [S] Assume that Rc > 0. Then the heat kernel p t (x , y) 
of the equation p t = Ap t satisfy the following estimate: 

Pt(x,y) > 



(4^)"/2 

In [9], Hamilton proved a matrix version of (II. ip for the heat equa- 
tion. The following is a matrix version of (11. ip . 

Theorem 1.6. Let p t be a solution of the equation p t = Ap t +(S7U, V ' pt). 
Assume that 

(1) the sectional curvature of the Riemannian manifold M is non- 
negative, 
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(2) the Ricci curvature is parallel, 

(3) V 2 (AU + ||Vt/| 2 ) < k 2 g for some non-negative constant k. 



Then 



V 2 log pt + X - V 2 U > - - coth(kt)g. 



Once again, if the underlying manifold is M. n and U(x) = — ||x| 2 , 
then 

V 2 ^AU + ^\VU\ 2 ^J = k 2 g, 

and 

V 2 \og Pt + ^V 2 U = ~coth(kt)g. 

Therefore, the inequalities in Theorem 11.61 are equalities in this case. 
By setting U = and letting k — > 0, we recover 

Theorem 1.7. j9] Let pt be a solution of the equation p t = Ap t . As- 
sume that 

(1) the sectional curvature of the Riemannian manifold M is non- 
negative, 

(2) the Ricci curvature is parallel. 

Then 

„ 2l 1 
V logp t > -—g. 

The above results can be proved using Bochner formula and the max- 
imum principle. However, instead of Bochner formula, we will prove a 
general result (Theorem 12.11 and I2.3j) using a moving frame argument 
as in [15]. In section 3 and 4, we show that the above generalization 
of Li-Yau estimate and its matrix analogue are simple consequences of 

Theorem O and El 

Aronzon-Benilan estimate is a differential Harnack inequality for the 
porous medium equation 

In section 5, we will prove a generalization of Aronzon-Benilan estimate 
using Theorem 12.11 and 12.31 Finally, we will prove sharp Laplace and 
Hessian type comparison theorems for the cost function (16. 1ft in section 
6. 
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2. Preliminaries 

In this section, we state and prove general results which will be used 
in the next few sections. For this, we will introduce some notations. Let 
M be a n- dimensional compact manifold without boundary equipped 
with a Riemannian metric denoted by (•,•} or g. The corresponding 
Riemann curvature tensor is denoted by Rm. Let F be a function 
on the space of all n x n matrices. We assume that F is invariant 
under orthogonal changes of variables (i.e. F(0 T AO) = F(A) for each 
orthogonal matrix O). For each linear map W : T X M — > T X M of the 
tangent space T X M at a point x, we set F(W) = F(W), where W is 
the matrix with ij-th entry equal to (W(vi),Vj) and {vi, ...,v n } is an 
orthonormal frame at x. This is well-defined since F is invariant under 
orthogonal changes of variables. Note that this condition is not needed 
or can be relaxed when the tangent bundle TM of M is parallelizable. 
For instance, when the manifold is the flat torus, this condition can be 
completely removed. Finally, if u, v, and w are tangent vectors, then 
u (g) v denotes the linear map defined by u <8> v(w) = (v, w) u. 

The following is a generalization of the Li-Yau estimate [TU] . 

Theorem 2.1. Assume that there is a non-negative function b t : M — >■ 
R, a time dependent vector field Y t on M , and a fibre-preserving bundle 
homomorphism W t : TM — > TM of the tangent bundle TM such that 

(1) F'(A)(B 2 ) > kiF(B) 2 for some non-negative constant k 1} 

(2) F'(VX)(W t + Rm(-,X)X) > k 3 for some constant k 3 , 

(3) F>(VX)(V(X t + V Xt X t ) + W t ) + k 2 F(VX t ) 2 
< F'(VX)(b t V 2 (F(VX t )) + V(F(VX t )) ® Y t ), 

(4) h + k 2 > 0, 

Then 

F(VX t ) < — - — a( kl+k2 ) k3 (t), 

where 

{VKcot(VKt) ifK>0 
\ ifK = 

^/^K coth(^/^K t) if K < 0. 

Remark 2.2. Note that the above theorem can be further generalized 
to include situation considered in [3] if F is allowed to depend on X, 
not just VX. However, we will not pursue this here. 

A matrix version of Li-Yau estimate was done by Hamilton [9]. The 
following is the corresponding matrix version of Theorem 12.11 
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Theorem 2.3. Assume that there is a non-negative function b t :M^- 
R, a time dependent vector field Y t on M , and a fibre-preserving bundle 
homomorphism W t : TM — y TM of the tangent bundle TM such that 

(1) w (X t ,w) is a closed 1-form, 

(2) Wt + Rm(-, X t )X t > k$g for some constant k%, 

(3) (v v (X t + VXt(Xt)), + k 2 (VX t (VX t (v)),v) + (W t v, v) 

< bt (AVX t (v),v) + (Vy t V w Xt, v) for each eigenvector of the 
linear map w \-y V w X t with the largest eigenvalue, 

(4) 1 + k 2 > 0, 

Then 

VXi < - —^- a {1+k2)k3 (t)g. 

As a consequence, we obtain the following estimate on the volume 
growth of a set under the flow of the vector field Xt if F — tr. 

Corollary 2.4. Under the assumptions of Theorem \2.1\ with F = tr, 

(b {kl+k2)k3 (t)y^vol(ipt(D)) 

is a decreasing function of time t, where 

(^sm(^Kt) ifK>0 
b K {t) = \t ifK = 

(^sinhi^Kt) ifK<0. 

The rest of this section is devoted to the proof of the above mentioned 
results. 

Proof of Theorem \2.1[ Let <p t be the one-parameter family of diffeo- 
morphisms defined by the vector field X t : (fit = X t (<*pt) and (fo{x) = x. 
Let 7(s) be a curve which start from x with initial velocity v: 7(0) = x 
and 7'(0) = v. Then 



Jt dipt{v) = T s Jt^ l{s)) 



s=0 



Let Vi(0), v n (0) be an orthonormal frame at a point x and let 
Vi(t), ...,v n (t) be the parallel transport of V\ (0) , . . . , v n (0) along the path 
(ft(x). Let A(t) be the matrix defined by 

n 
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It follows that 

n n 

{t)X t . 

i=i i=i 

Therefore, if ^-(t) = (V„ i(t) X t , then S(t) = A{t)~ l A{t) and we 

have 

S(f) = -A{t)^A(t)A(t)^A{t) + Aft)" 1 ^*) 

= -^(t) 2 + A(t)^ 1 i(t). 

On the other hand, if we differentiate the equation ip t = X t (tp t ), then 
we get 

^-<p t = X t {<pt) + V Xt X t (<pt) 



D D 

ds dt ' 



dt 

and 

^( 7 ( s ))| s=0 = + Vx * Xt 

By the definition of the Riemann curvature tensor Rm, it follows 
that 

D 2 / ■ \ 

jpd<p t (v) + Rm(dip t {v),X t {<p t ))X t {ip t ) = V dMv) \ X t + V Xt X t j . 

Therefore, by the definition of the matrix A(t), the following holds 
A(t) + A(t)(R(t) - M(t)) = 0, 

where 

Rij(t) = (Rm(^(f),X 4 (^(x)))Xt(^(x)),^(t)) 

and 

Mij(t) = (v Mt) [x t + Vx t -Xt) ,Vj(t) t 
By combining this with ( 12. II) . we obtain 

(2.2) 5(t) + 5(t) 2 + J R(t) = M(t). 
It follows that 

(2.3) j/(S(t)) + F'(S(t))(S(t) 2 + R(t)) = F'(S(t))(M(t)). 

Let t be the first time where F(VX to (<p to (x))) = /ca^(t ) for some 
point x, where k > 0. By assumption, we have 

(2.4) F'(VX t0 )(V(X t() + V Xt X t0 ) + W t ) + k 2 F{VX t0 ) 2 < 

at <pto(x). 

In the matrix notation, we have 

F'(S(t ))(M(t ) + W(t )) + k 2 F(S(t )) 2 < 0, 



s 



PAUL W.Y. LEE 



where W(t ) be the matrix with ij-th entry equal to (Wt(vi(t)),Vj(t)). 
By combining this with (12.31) and using the assumptions, we obtain 

j t F(S(t )) + hF(S(t )) 2 + k 3 

< j/(S(t )) + hF(S(t )) 2 + F'(S(t ))(R(t ) + W(to)) 

< ^-F(S(t )) + F'(S(t ) 2 ) + F'(5(t ))( J R(to) + mo)) 
at 

= F'(S(t ))(M(t ) + W(t )) 

< -k 2 F(S(t )) 2 . 

By the definition of to, we have kaxih) = F(S(to)) and kaxito) < 
j^F(S(t )). Therefore, the above inequality becomes 

ka K (t ) + (ki + k 2 )k 2 a K (t ) 2 + A; 3 < 0. 

Since satisfies 

(2.5) a K + a\ + K = 0, 

it follows that 

k((ki + k 2 )k - l)a(t ) 2 + k 3 - kK < 0. 

Therefore, we obtain a contradiction if k = kl + k2 and K < {k\ + 
k 2 )k 3 . Hence 

for all K < {ki + k 2 )k 3 . By letting If — >• (&1 + k 2 )k 3 , we obtain 
F(VX f ) < - — -^-a (fel+fe2 ) fe3 (t). 

M + K 2 

□ 



Proof of Theorem \2.3l Here, we use the same notations as in the proof 
of Theorem 12. 1 L By assumption the one-form v i— > (X t ,v) is closed. 
This is equivalent to (V v X t ,w) = (v,V w X t ). It follows that the ma- 
trices S(t) are all symmetric. Let to be the first time such that there 
is a point x and a unit tangent vector v in the tangent space T Vt MM 
at (ft(x) such that (V v X to ,v) = (S(t )v,v) = /ca^(to)- Here v denotes 
both the vector v and its matrix representation with respect to the 
orthonormal frame Vi(t), ...,v n (t). In particular, kaxito) is the largest 
eigenvalue of S(t ) with eigenvector v. By parallel translating along 
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geodesies, we extend v to a vector field still denoted by v. It follows 
that Vi> = and Av = 0. Therefore, the following holds by assumption 



In terms of the matrix notations, the above inequality becomes 



((M(t ) + k 2 S(t ) 2 + W(to))v, v) < 0. 
This, together with (12.2ft and ( 12.5ft . gives 

0<^-({S(t)v,v)-ka K (t)) 

at t=t 

= - (S(t ) 2 v, v) + ((M(t„) - i?(*o)K u) + A;a^ (to) 2 + kK 
< -(1 + fc 2 ) (S(t ) 2 v, v) - ((W(t ) + #(*o)K «> + fcM*o) 2 + fclf. 
By assumption, W(£) + R(t) > k 3 I. It follows that 



Therefore, we obtain a contradiction if k = and K < k 3 (l + k 2 ). 
It follows that 



(v,(X t0 + VX t0 (X t0 )),v) + A; 2 (VX (0 (VX i0 (i;)),t;) + (W to v,v) 

< b t0 (AVX t0 (v) lV } + (V Yto V v X to ,v) 

< b t0 A (VX t0 (v), v) + V YtQ (V v X t0 ,v) < 0. 



k(l - (1 + A; 2 )A;)ax(to) 2 + fcAT > fc 3 . 



1 + A; 2 



0- 



□ 



Proof of Corollary [O If F(VJC) > tr(VX), then 



-logdetA(t) <F(VX„ t(!E) ). 



It follows that 




where 





if X > 



if K = 



□ 
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3. A GENERALIZATION OF Ll-YAU ESTIMATE 

In this section, we prove the following generalization of the Li-Yau 
estimate. 

Theorem 3.1. Let p t be a positive solution of the equation 

(3.1) p t = Ap t + (Vpt, Wi> + U 2 pt. 

Assume that 

(1) Rc > 0, 

(2) A (-AC/i - H Wi| 2 + 2U 2 ) > k 3 . 
Then 

2Alogp t + AL^i > -nofca,(t). 

n 

Remark 3.2. If the underlying manifold is IR n , U\( I , and 

U 2 = 0, then (jl.3p satisfies the equality case 

-2Alog/9 t - At/i = nafcjt) 

of Theorem 13.11 

By integrating the above generalization of Li-Yau estimate, one ob- 
tains a Harnack inequality. For this, we need to consider the following 
functional 

ri|7(r)| 2 + y( 7 (r))dr, 

Js Z 

where 7 : [s , si] -> M and V = A£/i + ±| Wi| 2 - 2C/ 2 - 
Let c S0)Sl be the corresponding cost function defined by 

(3.2) c SOjS1 (x,y) = m{ H l| 7 (r)| 2 + V(j(r))dr, 

J s z 

where the infimum is taken over all paths 7 satisfying 7(^0) = x and 

t(*i) = y- 

Corollary 3.3. Under the assumptions of Theorem \3.l\ the following 
holds 

( ^ — — 

p S0 (x) V6*a(s ) / V 2 J 

If we let pt be the fundamental solution p t (x,y) of the equation 
(13.11) and let s — > in Corollary 13. 3[ then we obtain the following 
generalization of Cheeger-Yau estimate [5]. 
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Corollary 3.4. Let p t be the fundamental solution of the equation 
( Iff. 1]) . Under the assumptions of Theorem \3.1\ the following holds 

Pt(x,y)> (47rb^_(t)^ 2 exp ( -- (c , t (x,y) + Ut(y) - Ui(x)) 

Finally, we will show that the equality case in Corollary l3.4l is achieved 
by (jl.3p . More precisely, 

Theorem 3.5. Let p t be defined by A 1.3]) , U\(x) = — ||x| 2 , and U 2 = 0. 
Then 

p t (0,x) = exp (~ (cb, t (0,x) + U^x) - Z7i(0))J (47r6_ fc2 (t))-5 . 

Proof of Theorem \3.1[ If we specialize Theorem 12.11 to the case where 
F — tr and X t = Vh t , then the assumptions of Theorem 12.11 are 
satisfied if h = ±, tr(W t ) + TLc(X t , X t ) > k 3 , and 



A ^ + -|V/^ j + ^(A/^ + tr(Wi) 

<& t AA/i t + (VA/i t ,y t ). 
Let /i( = — 2 log p t — U\. Then the following holds 

ht + ^|V/i t | 2 = Aht + At/! + ^Vf/il 2 - 2C/ 2 . 

Therefore, under the assumptions of the theorem, (I3.3P holds with 
k 2 = and bt = 1. Hence, the result follows from Theorem 12.11 □ 

Proof of Corollaru \3.3\ Let 7 be a minimizer of (13. 2 p which satisfies 
7(so) — and 7(si) = £i- Using the notations in the proof of Theorem 
13. 1[ we have 

< J t Hl(t)) - (V^( 7 (t)),7(*)> + ^|Vfct|? (t , 

= Afc( 7 (f)) + AUMt)) + - 2U 2 (j(t)) 

<na^(t) + V(i(t)). 

n 

In the last inequality above, we have used Theorem 13.11 

By integrating the above inequality and noting that bx = bxCLK, we 
obtain 

/i Sl Oi) - h so (x ) < c S0|S1 (x ,xi) + nlog(6fca(si)) - nlog(6^(s )). 
By taking exponential of the above inequality, the result follows. □ 
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Proof of Corollary \3.4\ By Corollary 13.31 we have 

} > exp (-1 (c s>t (x,y) + UM - U x {x))^ 



Ps(x,X) V 2 J \ b t3\ s 

Since lim s ^o(47rs) n/ ' 2 p s (a;, x) = 1 (see [8]), the above inequality 
gives 

Pt(x,y) > exp (c ,t( x ^y) + ^lfa) - Uxix))^ (Anb^t) 

as claimed. □ 

Proof of Theorem \3.5l In this special case, the cost function (13. 2 p is 
given by 

(3.4) co, t (0,y) = inf J* ±\j(s)\ 2 + V(i{s))ds, 

where V(x) = —kn + |/c 2 |x| 2 and the infimum is taken over all paths 
7 satisfying 7(0) = and 7(t) = y. 

If x(-) is a minimizer of the above infimum, then it satisfies the 
following equations (see [7]) 

x = p, p — k 2 x. 

Since x(0) = and x(t) = y, it follows that 

. , sinhffcs) . , 
smh{kt) 

If we substitute this back into (13.4)) . then we obtain 

. . k\y\ 2 coth(H) 

coA Q >y) = 2 

A computation shows that 

p t (0,y) = exp ^-1 fcb, t (0,y) - ^M 2 )) (47r6_ fc3 

/ -k\y\ 2 \ /2yr(exp(2tfc) - 5 
~ 6XP \2(exp(2fci) - 1) J V A;exp(2H) J 

as claimed. □ 
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4. A GENERALIZATION OF HAMILTON'S MATRIX Ll-YAU ESTIMATE 

In this section, we show that the following generalization of Hamil- 
ton's matrix Li-Yau estimate is a consequence of Theorem 12.31 

Theorem 4.1. Let p t be a solution of the equation A3. Assume that 

(1) the sectional curvature of the Riemannian manifold M is non- 
negative, 

(2) the Ricci curvature is parallel, 

(3) -V 2 (A£7i + \\VU X \ 2 - 2U 2 ) > k 3 g, 

Then 

-2V 2 log A - V 2 fA < a k3 (t)g. 

Remark 4.2. If the underlying manifold is IR n , Ui(x) = — ||x| 2 , and 
U 2 = 0, then 

-V 2 (aCA + ilVf/il 2 - 2U^j = -k 2 g. 

If p t is given by (11 .4ft . then 

-2V 2 logp t - V 2 ^ = a k3 (t)g. 
Proof. We need the following lemma. 

Lemma 4.3. Assume that the sectional curvature of a Riemannian 
manifold is non-negative at a point x and the Ricci curvature Rc sat- 
isfies VRc x = 0. Then for any smooth function f , the following holds 

A(V v df(v))>(V v VAf,v). 

Here we consider Hessian of f as a self- adjoint operator on T X M and 
v is an eigenvector of this operator with the largest eigenvalue. 

Proof. Let v be a tangent vector at x and let e±, e n be an orthonormal 
frame at x. Let us extend v and to a vector field defined locally near 
x by parallel translation along geodesies. It follows that Vv(x) = and 
V ei V ei v (x) = (throughout this proof we sum over repeated indices 
without mentioning). Therefore, 

A(V v df(v)) = V ei V ei V v df{v). 

Let a be a 1-form and {3 be a two form. By Ricci identity, we have 

(1) V vl V V2 a(v 3 ) = V V2 V Vl a(v 3 ) - «(Rm(«i, v 2 )v 3 ), 

(2) \7 m V vl V V2 a(v 3 ) 

= V V4 V V2 V vl a(v 3 )-V VA a(Rm(vi,v 2 )v 3 )-a(V V4 Rm(vi,v 2 )v 3 ), 

(3) S7 vl V V2 P(v 3 ,V4) 

= V V2 V Vl P(v 3 , v A ) - /?(Rm(vi, v 2 )v 3 , v 4 ) - 0(v 3 , Rm(v 1 ,v 2 )v 4 ). 
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It follows that 

A(V v df(v)) = V ei (V v V ei df(v) - df(R m ( ei ,v)v)) 
= V ei V„V„d/( ei ) - V ei t^(Rm(e i ,v)v) 

- df(S/ ei Rm(ei,v)v) 

= V,V ei V^/( ei ) - V ei t^(Rm(e i ,v)v) 

- V v df(Rm(ei,v)ei) - V ej (J/(Rm(e j , v)v) 

- c//(V ei Rm(ei, v)v) 

= V v V v V ei df(ei) - V v df(Rm(ei,v)ei) 

- d/(V„Rm(e„ vfe) - V^df (Rm(e i; v)v) 

- V„(i/(Rm(e„ v)ei) - Ve j rf/(Rm(e i , v)v) 

- df(V ei ~Rm(ei,v)v) 

= (V„VA/,v) - 2V„rf/(Rm(e i ,t;)e ! ) 

- rf/(V„Rm(e„ u)e») - 2V ei rf/(Rm(ei, 

- df(V ei Bxn(e i ,v)v). 

Since the Ricci curvature is parallel, we have, by the contracted 
Bianchi identity, 

A{V v df{v)) = (V v VAf,v) - 2V„d/(Rm(e i ,v)e i ) - 2V ei df{Rxa{e h v)v). 

If ej is an eigenvector of the hessian of / with eigenvalue Aj and v is 
an eigenvector of the hessian of / with the largest eigenvalue A, then 

A(V v df(v)) 

= (V„VA/, v) + 2ARc(t;, v) - 2\ (a, Rm(ei, v)v) 
> (V v VAf,v). 

Here we use the assumption that the sectional curvature is non-negative. 

□ 

When X t = Vh t , the conditions become W(t) + R(t) > k 3 I and 

V„V (ht + \\Vh t \ 2 ^ ,v\+k 2 ((V 2 h t ) 2 v, v) + (W t v, v) 

< b t (AV 2 h t (v),v) + (V Yt V 2 h t (v),v) 

for each eigenvector v of the symmetric operator V 2 h t with the largest 
eigenvalue. 

Recall that if p t is a positive solution of the equation 
Pt = Ap t + (Vpt, VC7i) + U 2 p t , 
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then h t = —2 log p t — U\ satisfies 

ht + \\Vh t \ 2 = Ah t + Af/x + ^Vf/il 2 - 2U 2 . 
It follows that 

V 2 (h t+ l -\Vh t \ 2 \ +W t = V 2 Ah t , 

where W t = -V 2 (AC/i + ||Wi| 2 - 2t/ 2 ) . 

Therefore, if we assume that the Ricci curvature is parallel, the sec- 
tional curvature is non-negative, and Wt > k^g, then 

V 2 (ht + ^|V/it| 2 ^ (v),vj + (W t (v),v) < (AV 2 h t (v),v) . 

It follows that 

V 2 h t = — 2V 2 log p t - V 2 ^ < a k3 (t)g. 

□ 



5. A GENERALIZATION OF ARONZON-BENILAN ESTIMATE 

Aronzon-Benilan estimate [2] is a differential Harnack inequality for 
the porous medium equation 

In this section, we apply Theorem 12.11 and prove the following gen- 
eralization of Aronzon-Benilan estimate. 

Theorem 5.1. Let p t be a smooth positive solution of the equation 

Pt = A(pD + u P 2 ~ m . 

Assume that Rc > 0, AU > andm-l + ^>0. Then 

— 7 — 2 m 7 n 

^2/fY\j . , — 2?^ , , 

-A{p t ) < — — q fc 3 (2 + n(m-l)) (t). 

m — 1 2 + n(m — 1) 2„ 

Proof. A computation shows that h t = yz^pT' 1 satisfies 

ht + \\^h t \ 2 = ^(1 - m)h t Ah t - 2mU. 

It follows that 

A (ht + ^\Vtk\ 2 j + l -[m - \){Ah t ) 2 + 2mAU 

= (1 - m) (Vh t , VAh t ) + mp™- 1 A Aht. 
The rest follows from the assumptions and Theorem 12.11 □ 
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6. On Laplacian and Hessian comparison type theorems 

In this section, we prove versions of Laplacian and Hessian type 
comparison theorems for the following cost function 

(6.1) c s , t (x,y)= inf f \\^r) - W^r))] 2 - U^{r))dr. 

7(s)=z,7W=2/ J s 1 

More precisely, 

Theorem 6.1. Assume that 

(1) the Ricci curvature of the underlying manifold M is non-negative, 

(2) A (U2 — llVL^ij 2 ) > k% for some negative constant k%. 
Then the cost function co,t defined by $1.60 satisfies 



A x c , t (x ,x) < ^-k 3 ncoth ( \/-—t 



n 



wherever c 0i t(x , •) is twice differentiable. 

Theorem 6.2. Assume that 

(1) the sectional curvature of the underlying manifold M is non- 
negative, 

(2) V 2 (U2 — llVL^il 2 ) > k 3 g for some negative constant k 3 . 
Then the cost function co,t defined by $1.6\) satisfies 

V 2 c 0i t(x ,x) < v^acoth (y/^fat\ g 

wherever Coj(%o, ') twice differentiable. 

Remark 6.3. The function x t— > coj(%o,%) is locally semi-concave. In 
particular, it is twice differentiable Lebesgue almost everywhere by 
Alexandrov's theorem. Therefore, the conclusions in Theorem 16. II and 
16.21 hold Lebesgue almost everywhere (see [15] for the definitions and 
the results). 

Remark 6.4. We can see that the above theorems are sharp by looking 
at the case M = M™, U Y = 0, and U 2 (x) = -^\x\ 2 . We have 

V 2 (u 2 - i|Wi| 2 l = ~k 2 I, A (u 2 - -|V^i| 2 l = ~k 2 n 



which is the equality case in the second conditions of Theorem 16.11 and 
6.2[ We also have 

V 2 c 0i t(0,x) = kcoth(kt)I, A x c j(0,x) = kncoth(kt) 

which gives the equality case in the conclusions of Theorem 16.11 and 
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Remark 6.5. A Bishop-Gromov type volume comparison theorem fol- 
lows from Corollary 12.41 

The proof of Theorem 16.21 is similar to that of Theorem 16.11 and will 
be omitted. 

Proof of Theorem \6.1[ If Cqj{xq,x) is smooth, then the result follows 
from Theorem 12.11 and 12. 31 Indeed, the Legendre transform of the 
Lagrangian 

L(x, v) = l - gij (x)(v l - g il {x){U x ) Xl {x)){v^ - g' k {x)U Xk {x)) - U 2 (x) 

is given by 

H(x,p) = sup [p{v) — L(x, v)] 

v&T x M 

= ^9 ij {x)PiPj + 9 ij (x)Pi(Ui) Xj (x) + U 2 (x). 

Here we sum over repeated indices. 

The corresponding Hamilton- Jacobi equation is given by 

(6.2) f t + ^\Vf t \ 2 + (VU 1} V/ t ) + U 2 = 

and c 0jt (x ,x) is a particular solution (see [1]). 
If we set X t = V (co,t(xo, •) + U\), then 

Vx t X t ) - i|V 



2' 



Therefore, 



A x c 0it (x ,x) < naks(t) = \/ -k z n coth [ \/-—t 



n 



by Theorem 12. 1[ 

In general, if X IS cL point where c Qjt (x , ■) is twice differentiate, then 
there is a unique minimizer 7 to the infimum (jl.6p joining Xo and x. 
Moreover, c 0)S (x ,-) is smooth at 7(3) for each s in (0,t) (see [1]). 
Therefore, the proof of Theorem 12.11 still applies. Note that, in this 
case, (12.41) is an equality. □ 
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